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Abstract 

Wigner crystallization in mesoscopic quantum dots containing only few (N < 50) 
electrons exhibits a number of interesting peculiarities: (i) there exist two distinct 
crystal phases, and (ii) the phase boundary sensitively depends on the precise parti- 
cle number. In this paper we demonstrate that this behavior can be used to control 
the collective transport properties by adding or removing a single electron. 
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One of the most exciting properties of mesoscopic systems is the dependence of 
their properties on the particle number. Besides the familiar electron addition 
spectra in quantum dots, e.g. ([!]), recently, another example has been found, 
cf. (0; and references therein: Wigner crystallization in mesoscopic electron 
clusters in two-dimensional quantum dots. Simulations of a small number of 
electrons confined in a spherically symmetric harmonic trap and interacting via 
the Coulomb potential £|) revealed that the location of the phase boundary 
of the Wigner crystal in the density-temperature plane sensitively depends on 
the precise particle number. In this paper we discuss principle possibilities of 
taking advantage of this behavior for applications. 

The particle number sensitivity originates from the configuration symmetry 
of the cluster ground state. Minimization of the total energy computed from 
the hamiltonian 



i=l * m i i=l z i<j 
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(m* and e& are the effective electron mass and background dielectric constant, 
respectively), yields a spherical shell structure details of which vary strongly 
with the particle number, see Fig. 1. "Magic" clusters (those with an integer 
ratio of particle numbers on the outer and inner shell - e.g. N = 19) show 
a clear hexagonal crystal structure (Fig. 1, left part). In contrast, non-magic 
ones (e.g. N = 20) are dominated by the spherical trap symmetry. This has a 
strong effect on the stability of the crystal phases, see Fig. 1, right part. 
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Fig. 1. Left two figures: Ground state configuration of the "magic" cluster N = 19 
(top) and N = 20 (bottom) in the cluster plane. Each dot is an electron. 
Right figure: Phase diagram of the Wigner crystal of 2D clusters with N = 19 and 
N = 20 electrons. The outer (inner) lines are the radial (angular) melting phase 
boundaries. The dotted diagonal line separates the classical (left) from the quantum 
(right) crystal. The classical part of the boundaries is given by lines with constant 
r, with the values (from top to bottom) Tffi = 83, Tf 9 M = 154, T? g M = 330 and 



rg/* = 3.4-10 u . 
and T = 1/r. 



The dimensionless density and temperature are given by n 
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Generally, crystal-like behavior (^) is found below a critical temperature (of 
the order of a few K in semiconductor quantum dots) in a finite density 
interval, see Fig. 1. If the density is reduced below a critical value ni(T, N), the 
system undergoes a transition to a state resembling a classical liquid. Similarly, 
above a second critical density ri2(T,N) (with ri2 > n\) a Fermi liquid-like 
state (in mesoscopic systems, sometimes referred to as Wigner molecule) is 
reached. While this general behavior is analogous to macroscopic systems, 
e.g. (|jp, crystallization in few-electron systems shows a number of interesting 
peculiarities (|7|): (i) strong N dependence of the phase boundary and (ii), 
the existence of a second phase boundary which is embedded into the first 
one where the crystal structure is transformed from a completely ordered 
( "OO" , fully localized electrons) state into a partially ( "RO" , radially) ordered 
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one where electrons are rigidly confined to crystal shells which, as a whole, 
can rotate against each other. The phase boundary is determined by critical 
values of the coupling parameters: the low-density boundary by the classical 
parameter, T = e 2 /(ebr k B T), and the high-density limit by r s = r /a,B, cf. 
Fig. 1. Here, as is the effective Bohr radius, clb = h 2 eb/(me 2 ), and r is 
the mean interparticle distance, approximately given by the balance of the 
repulsive Coulomb force and the radial confinement force of the harmonic 
potential with strength u : e 2 /(e?> r o) = rnu} 2 r . A particularly strong N 
dependence is observed for the "00" phase: "magic" clusters (N = 12, 19 
etc.) are found to be much more stable than clusters having one electron more 
or less (@; |]). For example, the orientational melting parameters T OM and 
Tg of the cluster with N = 19 are 9 orders of magnitude lower than those 
for N = 20, see Fig. 1. 

This behavior can be exploited for a non-traditional control of crystallization. 
In addition to changing temperature or/and density (confinement strength), 
crystallization in mesoscopic systems can be achieved by variation of the par- 
ticle number without change of T and n. For example, choosing a point in 
the temperature-density plane, Fig. 1, which is located between the radial 
(orientational) melting curves of N = 19 and N = 20 and switching between 
the two particle numbers is equivalent to a crossover between crystal-like and 
liquid-like (00 and RO crystal) behavior. We will demonstrate this below for 
the orientational melting curves in the classical part of the phase diagram, i.e. 
we fix the classical coupling parameter T between rf? M and rf 9 M , (see Fig. 1, 
right part). The main difference between the orient at ionally ordered (with 
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Fig. 2. Spontaneously excited angular currents of the inner ("in") and outer ("out") 
shell of an N = 20 cluster at T = 500 -C T^q 1 averaged over To = 2it/ujq. The outer 
shell is pinned. The thick dashed line is j ou t averaged over 4Tq. 

N = 19) and disordered (N = 20) state will be that the latter is able to 
support inter-shell rotational excitations. This should give rise to macroscopic 
currents and magnetic fields. To estimate this effect, we compute the total 
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angular current created by N electrons confined to one of M shells of radius 
Rk rotating with angular frequency uj k (positive or negative depending on the 
direction), 1$ = eY^kLi N k uJkRk- The associated magnetic field on the sym- 
metry axis in a distance z above the cluster plane is directed normal to the 
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plane: B z (z, t) = f Eifeli N k u k (t) [1 + (z/R k ) \ For example, for N = 20, 
we have N x = 7, N 2 = 12, R 2 ^2R 1 ^2- r , so that If « er (7 • u x + 24 • u 2 ). 
The relation between uj\ and u 2 depends on the excitation conditions. For 
example, if the total angular momentum of the excitation is zero, one readily 
finds uj\/uj 2 —4N 2 /Ni, i.e. the inner shell rotates approximately 7 times 
faster. Alternatively, rotation of some shells may be inhibited due to defects 
(pinning). For definiteness, in the following we consider the cluster A^ = 20 
with the outer shell pinned. 

At a given coupling T and confinement energy muJorl/2, thermal fluctuations 
spontaneously excite rotational and vibrational degrees of freedom both, in 
radial and angular direction. The thermal energy per particle is ksT and is 
T times smaller than the confinement energy. As our simulations show (see 
below), approximately half of it is converted into angular kinetic energy (ro- 
tations and oscillations), i.e. E rot m \muj 2 R\/T . Depending on the excitation 
conditions, this gives rise to rotations of the inner shell electrons with U\ max ~ 
luq/V2T. Inserting this result into the expression for B z and averaging over 

-3/2 



fie 2 Ni 



,2 



the fast vibrations, we obtain {B z (z)) max ~ , 1 + (z/roY 

2e b r ^mr T 1 

Interestingly, the current and the magnetic field increase as the square root of 
temperature. In contrast, the cluster of 19 electrons does not support intershell 
rotations. Here, the excitation energy is completely converted into vibrations. 

We have verified this concept by performing classical molecular dynamics sim- 
ulations for N = 19 and A^ = 20 in a wide range of T values (§ [|). As expected, 
the cluster with A^ = 19 (pinned to suppress the trivial rotation of the whole 
system) supports, at T — 500 > rf 9 M , no intershell rotations, and only vibra- 
tions are excited, even for purely rotational initial fluctuations (0). In contrast, 
for A^ = 20, strong collective rotational motion of the inner shell electrons is 
excited, see Fig. 2. In all cases, we found that, in the spectrum, the vibra- 
tional excitations are well separated from the rotations - the former occur at 
significantly higher frequency, dominantly at frequency u (and uj /A). After 
averaging over 4T , with T = 2tt/u , a "persistent" collective current (jj n ) of 
the inner shell is observed which exceeds the averaged fluctuating signal (j ou t) 
of the pinned outer shell by at least two orders of magnitude. The direction of 
the shell rotation and its slow time-dependence of (j in ) are determined by the 
(random) excitation conditions and is not relevant. As an independent test we 
computed the potential barriers for intershell rotation using Quantum Monte 
Carlo simulations (^). These simulations effectively average over the random 
excitations and yield, for A^ = 20 at T = 500, potential barriers which are 
practically zero. 
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For practical applications, one may think of designing a suitable external ro- 
tational excitation in such a way that it is sufficiently weak so it does not 
overcome the barrier of the N = 19 cluster. If applied to the N = 20 clus- 
ter, the same excitation will easily give rise to intershell rotations. Such an 
excitation could be e.g. a constant or pulsed circularly polarized electric field. 
Based on our simulations, we expect that the response of a (pinned) N = 20 
cluster will be a nearly dissipationless rotation of the inner shell which stops 
after removal of one electron @. This allows for a completely new kind of 
single-electron controlled devices. Such collective rotations of a group of elec- 
trons can be easily transmitted, e.g. in a system of multiple layers ( |T0"D and 
may find applications for quantum computing. Naturally, realization of this 
concept will require very clean samples, and detection of the weak circular 
currents is another important issue to be solved. At the same time, high sen- 
sitivity of these few-electron clusters to very weak rotational excitations could 
be of interest for applications by itself. 
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